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ABSTRACT 


In this paper, we investigated the issue of black hole masses and minimum 
timescales of jet emission for blazars. We proposed a sophisticated model that 
sets an upper limit to the central black hole masses M, with the minimum 
timescales of variations observed in blazars. The value of presents 

an upper limit to the size of blob in jet. The blob is assumed to be generated in 
the jet-production region in the vicinity of black hole, and then the expanding 
blob travels outward along the jet. We applied the model to 32 blazars, 29 of 
which were detected in gamma rays by satellites, and these are on the order 

of hours with large variability amplitudes. In general, these M, estimated with 
this method are not inconsistent with those masses reported in the literatures. 
This model is natural to connect M, with Atj*^ for blazars, and seems to be 
applicable to constrain M, in the central engines of blazars. 

Subject headings: black hole physics - BL Lacertae objects: general ~ galaxies; active 
- galaxies: jets - quasars; general 
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INTRODUCTION 


Blazars are radio-loud active galactic nuclei (AGNs), including BL Lacertae objects 
(BL Lac objects) and flat spectrum radio quasars (FSRQs), characterized with some special 
observational features, such as luminous nonthermal continuum emission from radio up 


to GeV/TeV energies, rapid varia 
of their compact radio cores (e.g. 


jility with large amplitu des, and superluminal motion 


Urrv fc Padovaui 


199511 ■ These unusual characteristics 


originate fro m the Doppler boosted em ission of a relativistic jet with a small angle to the 


line of sight flBlandford fc Rees 


19781) . The intranight or intraday variability (IDV) is an 


i ntrinsic phenomenon, an d tightly constrain the diameters of the emitting regions in blazars 
flWagner fc Witzel Ill995h . The emitting region in the relativistic jet is usually simplified as 
a blob. The size D of blob can be limited hj D < hAtobc/(l -|- z), where 6 is the Doppler 
factor of jet, Atob is the observed timescale of variations, is the redshift of source, and c 
is the speed of light. The IDV limit D to be smaller than the size of solar system by At < 
1 day in the source rest frame. These timescales of the variations with large amplitudes in 
the optical-gamma-ray bands might have some nnderlying connection with the black hole 
masses of the central engines in blazars. 

The relativistic jets can be generated from inner accretion disk in the vic i nity of black 
hole (e.g. 


Penrose 

1969; 

Blandford & Znaiek 

1977; 

Blandford & Pavne 

1982; 

Meier et ah 


200ll). Observations show that dips in the X-ray emission, generated in the central 


and microq uasars (e.g. 


20091, 


elections of bright superluminal radio 

mots in the lets of AGh 

Marscher et ah 

2002; 

Arshakian et ah 

2010; 

Ghatteriee et ah 


201ll ). The dips in the X-ray emission are well correlated with the ejec tions of 


brig ht superluminal knots in the radio jets of 3G 120 (IGhatterjee et ah 


111 flGhatteriee et ah 


20091) and 3G 


201ll) . An instability in accretion flow may cause a section of the 


inner disk break off, and the loss of this section leads to a decrease in the soft X-ray flux. 
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observed as a dip in the X-ray emission. A fraction of the section is accreted into the 
event horizon of the centra black hole. A considerable portion of the section is ejected into 
the jet, observed as the appearance of a superluminal bright knot. General relativistic 
magnetohydrodynamic simulations showed production of relativistic jets in the vicinity of 
black hole with jet-production region around 7-8rg for the Schwarzchild black hole and 
being of the order of the radius of the ergo sphere = 2r„ for the Kerr black hole, where 
Tg is the gravitational radius of black hole flMeier et ah II 2 OOII) . For the Schwarzchild black 
hole, the inner radius of the accretion disk is about equal to that of the marginally stable 
orbit 6rg, which is comparable to the size of the jet-production region. Thus the initial 
size of a blob emerged from the jet-production region may be comparable to that of the 
marginally stable orbit or that of the ergosphere. The blob will expand as it travels outward 
along the jet. When the blob pass the site of dissipation region in the jet, it will produce 
the corresponding variations in the optical-gamma-ray regimes. The minimum timescales 
of the variations are likely to be related with the masses of the central black holes. Then 
the black hole masses M, can be constrained with the observed minimum timescales 
of variations in the optical-gamma-ray regimes. In this paper, we attempt to construct a 
model that constrains M, with for blazars. 


The structure of this paper is as follows. Section 2 presents method. Section 3 presents 
applications. Section 4 is for discussion and conclusions. 


2. METHOD 


Assuming the blob has a size of Dq in the jet-production region (with a size comparable 
to Do), and a size of Dr at the location Djet in the jet from the central engine (see Figure 
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1), we have an equation between Dq and Dr 

Dr = Do + 2nexp^ = Do + (1) 

^jet Wjet 

where Uexp is the average expansion velocity of blob in the jet between the central engine 
and the location Djet, and Ujet is the corresponding average bulk velocity of the blob. The 
jet velocity fJjet is in a relativistic region. The expansion velocity Texp is not in a relativistic 
region. Then Texp -C Tjet, and we have 

Oo < Cr. (2) 

The size Dr of a blob at the location Djet can be constrained by the observed minimum 



Fig. 1.— Sketch of axial cross section of geometry. Dq is the size of blob in the jet formation 
region around the central engine, and Dr is the size of blob at the position Djet in the jet 
from the central engine. 


timescale of the variations from the blob, and we have 


Dr < 


1 + ^ 


(3) 
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where 6 is the Doppler factor, is the redshift of source, and c is the light speed. Combining 
equations (2) and (3), we have 


]j < 

1 + Z 


c. 


(4) 


The inner radius of accretion disk is usually taken to be around the marginally stable 
orbit of disk surrounding the central black hole. The radius of marginally stable orbit of 
disk is 

= ''"g{3 + -^2 T [(3 — ^i)(3 + Zi + 2 ^ 2 )]^'^^}, (5) 

where rg = GM,/(? is the gravitational radius of black hole with a mass of M, 


Zi = 1 + (1 — + (1 — and Z 2 = (3j^ + 1 Bardeen et ah 1972). 

Here, j = J/Jmax is the dimensionless spin parameter of black hole with the maximum 
possible angular momentum Jmax = GMl/c with G being the gravitational constant. In 
the case of the prograde rotation, r^s = 6rg for j = 0 and r^s = 'Tg for j = 1. General 


relativistic magnetohydrodynamic simulations sh ow that in t 


le Sc 


jet-production region has a size around 14-16rg flMeier et ah 


iwarzchild case the 


20011) . comparable to the 


diameter of the marginally stable orbit, = 12rg. In the Kerr case, the jet-production 
region must be of the order of the diameter of the ergosphere Dp = 2re = 4rg, where rg is 


the equatorial boundary of the ergosphere (IMeier et ah 


20011) . Thus the jet-production 


region span about from to for 0 < j < 1, and then we take Dq =4-12rg. From 
equation (4) and Dq =4-12rg, we have 


M. < 5.086 X IO^-^Mq (Do = Dg, j ~ 1), 

1 “r Z 

M. < 1.695 X 10"—(Do = D^g, j = 0), 

1 + z 

where is in units of seconds. Equations (6a) and (6b) can be unihed as 

M, < 1.695-5.086 x 10^- 


(6a) 

(6b) 


(7) 


1 + ^ 
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3. APPLICATIONS 


The new method is applied to 32 blazars. These blazars have variability timescales 
of the order of hours with large variability amplitudes and a redshift range from 0.031 to 
1.813. These timescales were reported in the literatures for the optical-X-ray-gamma-ray 
bands. The satellites detected GeV gamma rays from 29 out of 32 blazars. The details 
of these blazars are presented in Table 1. The minimum timescales of gamma rays are 
taken to be the doubling times of fluxes. The optical minimum timescales were reported 
in the literatures presented in column (6) of Table 1. The doubling times of X-ray fluxes 
are taken as the X-ray minimum timescales. These observed minimum timescales are 
listed in column (4), and their corresponding references are presented in column (6). The 
optical~ o^-rav emission are mo stly the Doppler boosted emission of jets for gamma-ray 


blazars flGliisellini et ah 


(iGhisellini et ah 


19981) . A value of 5 ~ 10 was adopted for GeV gamma-ray blazars 


2010 1. We will take 5 = 10 to estimate M, with formula (7). 


The Schwarzchild and Kerr black holes are considered in estimates of M, from formula 

(7) . The estimated black hole masses are denoted by and in columns (7) and 

(8) of Table 1, respectively. In the Schwarzchild case, spans from 10®'^^ to 10®-®® Mq. 

In the Kerr case, spans from 10®'®^ to 10^®-^® Mq. We compared these estimated M®J^® 

to those masses Mbh obtained with other methods reported in the literatures. In Figure 2 
of Mbh versus there are hve blazars below the line of Mvar = Mbh- Because M^jJ* is 

only an upper limit of M,, the area below this line is not allowed for these blazars. When 
the central black holes are the Kerr ones, these upper limits of masses are increased by a 
factor 3, and four out of hve blazars are moved into the area above the line of Mvar = VTbh- 
Only 4G -1-38.41 is just below the line of Mvar = Mbr. The populations of M, in Figure 

2 show that these Mvar estimated from formula (7) are reasonable upper limits on M, for 
blazars. 














Fig. 2.— Black hole mass Mvar versus Mbh estimated in other methods in the literatures. 
The solid line is M^ar = TTbh- The triangles are M^ar in the Schwarzchild case, and the solid 


circles are iW^ar in the Kerr case. 
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DISCUSSION AND CONCLUSIONS 


Morini et al. 


fll986l) thought X-ray emi ssion to be produc e d ver y close to the inner 


Aharonian et ah 


( 20071) limited the Doppler 


engine in BL Lacertae object PKS 2155-304. 
factor by the black hol e mass and th e varia bility timescale of very high energy gamma-ray 
flare of PKS 2155-304. 


Miller et al. 


(119891 ) reported the rapid variations on timescale as 


short as 1.5 hours for BL Lacertae in the optical flux, and the minimum timescale for the 
variations was used to place constraints on the size of the emitting region. They assumed 
that these variations were produced in the vicinity of a supermassive black hole, and then 
determined a black hole mass with the minimum timescale from a formula 


6G'(1 + z) ’ 

for the Schwarzchild black hole. At that time, they thought that relativistic beaming 


( 8 ) 


Xie et al. 

(2nn2c 

) deduced a formula from 

Abramowicz & Nobili 

- ^ 

00 


A+ob 

M. < 1.62 X 10^— 


1 + z 


Qi 


(9) 


which gives an upper limit on M,. These above two equations are based on assumption 
of accretion disk surrounding a supermassive black hole, and the optical flux variations 
are from the accretion disk. Obviously, these two formulae are applicable to estimate M, 
for non-blazar-like AGNs or some AGNs with weaker blazar emission component in fluxes 
relative to accreti on disk emission component. Gonsidering the relativistic beaming effect. 


Xie et al. 


(l2002d) deduced a new formula from formula (9) 


M, < 1.62 X 


1 + x 


0 - 


( 10 ) 


Formu la (10) was applied to blazars, especially BL Lac objects (see 


Xie et al. 


2002d, 


2005bl) . In this paper, we proposed a sophisticated model to constrain the black hole masses 
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using the rapid variations with large amplitudes for blazars. The model is suitable to 
constrain M, in blazars using the minimum timescales of variations of the beamed emission 
from the relativistic jet. Formula (7) is the same as formula (10) in except of the coefficients 
in the right of formulae, but these two corresponding models are essentially different in 
their origins. 


Spectral energy distributions of blazars consist of two broad peaks (e.g. 


Ghisellini et ah 


199811 . The hrst peak is produced by the synchrotron radiation processes of relativistic 


electrons in a relativistic jet. The second one is generally believed to come from the 
inverse Compton scattering processes of the same population of electrons. In general, the 
optical-gamma-ray emission are believed and/or assumed to be produced in the same 
region, simplihed as a blob in simulating the spectral energy distributions of blazars. The 
coincidence of a gamma-ray flare with a dramatic change of optical polarization angle 


provides evidence 


flAbdo et al. 


or co-spatiality of optical and gamma-ray emission regions in 3C 279 


201011 . Thus the timescales of the optical-X-ray-gamma-ray variations are 


used in formula (7). For this model, the blob size is assumed to increase linearly as it 
travels outward along the jet (see equation (1)). This linear growth assumption is only an 
approximation to the actual growth. The assumption will not change formula (2), and then 
don’t alter formula (7). 


The errors associated with this method are quite high. Besides the poorly known 
details of the model itself, related to our ignorance of the exact distance from the black 
hole where the emission is produced, a large uncertainty is related to the value of the 
adopted Doppler factor. A value of 5 = 10 is assumed, which may be in excess or short of 
the real value by at least a factor 3, resulting in a total uncertainty of at least an order of 
magnitude. In fact, it is possible for 32 blazars listed in Table 1. The Doppler factor 6 can 
be estimated from the Lorentz factor and the viewing angle of jet adopted to model spectral 
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energy distributions of bright Fermi blazars in 


Ghisellini et al. 


(j2O10). There are 21 Fermi 


blazars with new estimated 6, and these new 6 are from 13 to 28 with an average of 17. So, 
the upper limits of black hole masses are increased by a factor 1.3-2.8 when adopting these 
new 6 for the 21 blazars. Thus the 29 Fermi blazars out of 32 blazars we employed will 
have a similar case for the upper limits of black hole masses. The adopted Doppler factor 
of 5 = 10 may result in a large uncertainty with a factor 1.3-2.8 in the upper limits of M,. 
Another larger uncertainty with a factor 3.0 arises from the ignorance of spins of the central 
black holes in blazars (see formulae (6) and (7)). Thus the two large uncertainties will lead 
to errors of 3.9-8.4 in the upper limits of M, estimated with this method for blazars. 


In this paper, we proposed a sophisticated model to constrain the central black hole 
masses M, with the observed minimum timescales of variations in blazars. The size 

of a blob in the relativistic jet can be constrained with The blob is assumed to be 

ejected from the jet-production region in the vicinity of black hole, and then to expand 
linearly as it travels outward along the jet. The model is applied to 32 blazars, out of which 
29 blazars were detected in the GeV gamma-ray regime with the satellites. Their observed 
minimum timescales are on the order of hours with large variability amplitudes in the 
optical-X-ray-gamma-ray bands. In general, these M, estimated from by formula (7) 

are not inconsistent with those masses reported in the literatures. This indicates that this 
model is applicable to constrain M, in the central engines of blazars. This model is more 
natural to connect M, with for blazars. Due to the ignorance of the black hole spins 

and the real Doppler factors, the uncertainty of the upper limits of masses will be about 
3.9-8.4 in this method. The upper limits of masses will have uncertainty of 3.0 because the 
black hole spins are unknown for blazars, though we have the real Doppler factors rather 
than the adopted value of ^ = 10 in this paper. 
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Table 1. Sample of blazars 


Blazar name 

(1) 

z 

(2) 

Satelli.'* 

(3) 

A/ ■ ^ 

(4) 

Band 

(5) 

Ref. 

(6) 

mK-- 

(7) 

(8) 

Mbh" 

(9) 

Ref. 

(10) 

3C 66A 

0.444 

EG,LAT 

3.73 

V 

19 

9.28 

8.80 

8.60 

24 

AO 0235+164 

0.940 

EG,LAT 

4.74 

R 

9 

10.16 

9.68 

9.30 

24 

RBS 0413 

0.190 

LAT 

3.19 

V 

18 

8.82 

8.34 

7.95 

26 

PKS 0420-01 

0.916 

EG,LAT 

4.65 

o 

1 

10.07 

9.60 

9.03 

26 

PKS 0537-441 

0.894 

EG,LAT 

4.66 

o 

1 

10.09 

9.61 

9.20 

24 

PKS 0548-322 

0.069 

None 

3.14 

R 

21 

8.82 

8.34 

8.15 

22 

S5 0716+71 

0.300 

EG,LAT 

3.17 

V 

7 

8.76 

8.28 

8.10 

24 

PKS 0735+17 

0.424 

EG,LAT 

3.32 

B 

19 

8.87 

8.40 

8.40 

24 

OJ 248 

0.941 

EG,LAT 

3.65 

R 

20 

9.07 

8.59 

8.41 

27 

OJ 287 

0.306 

EG,LAT 

3.59 

R 

19 

9.18 

8.70 

8.51 

23 

MrK 421 

0.031 

EG,LAT 

3.26 

TeV 

4 

8.95 

8.48 

8.29 

22 

4C +29.45 

0.725 

EG,LAT 

4.40 

O 

17 

9.87 

9.39 

8.56 

27 

W Com 

0.102 

LAT 

3.58 

V 

16 

9.24 

8.77 

7.40 

24 

4C+21.35 

0.432 

EG,LAT 

3.92 

GeV 

3 

9.47 

8.99 

8.20 

24 

3C 273 

0.158 

EG,LAT 

4.33 

V 

11 

9.97 

9.50 

9.39 

25 

3C 279 

0.536 

EG,LAT 

3.38 

V 

19 

8.90 

8.42 

8.43 

26 

PKS 1406-076 

1.494 

EG,LAT 

4.76 

R 

14 

10.07 

9.59 

9.40 

24 

PKS 1510-089 

0.360 

EG,LAT 

3.39 

R 

20 

8.96 

8.49 

8.31 

27 

AP Lib 

0.049 

LAT 

2.95 

O 

1 

8.64 

8.16 

8.09 

22 

PKS 1622-29 

0.815 

EG,LAT 

4.14 

GeV 

8 

9.59 

9.11 

9.10 

24 

4C +38.41 

1.813 

EG,LAT 

4.76 

GeV 

2 

10.02 

9.54 

10.10 

24 

3C 345 

0.593 

None 

3.60 

I 

19 

9.10 

8.63 

8.43 

27 

MrK 501 

0.034 

LAT 

3.80 

R 

19 

9.49 

9.01 

9.20 

22 

3C 371 

0.051 

LAT 

3.38 

X 

12 

9.06 

8.59 

8.50 

22 

PKS 2005-489 

0.071 

EG,LAT 

3.84 

X 

5 

9.52 

9.04 

9.00 

26 

OX 169 

0.213 

LAT 

3.69 

O 

1 

9.31 

8.84 

8.14 

27 

PKS 2155-304 

0.116 

EG,LAT 

2.95 

O 

10 

8.61 

8.13 

7.10 

24 

BL Lacertae 

0.069 

EG,LAT 

3.30 

B 

15 

8.98 

8.50 

8.23 

26 

CTA 102 

1.037 

EG,LAT 

4.63 

O 

1 

10.03 

9.55 

9.10 

24 

3C 454.3 

0.859 

EG,LAT 

3.97 

O 

13 

9.41 

8.93 

9.17 

26 

OY +091 

0.190 

None 

3.30 

B 

15 

8.93 

8.45 

8.62 

26 






Table 1—Continued 


Blazar name 

z 

Satelli.^ 

At ■ 

Band 

Ref. 

MKe- 

MS-hc 

Mbh" 

Ref. 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

lES 2344+514 

0.044 

LAT 

3.70 

X 

6 

9.39 

8.91 

8.80 

26 


Note. — Column 1: Blazar names; Column 2: redshifts of objects; Column 3: The 7 -ray 
satellites detected 7 rays from objects; Column 4: the minimum timescales of variations of 
objects; Column 5: the bands where the minimum timescales measured, B: the B band; V: 
the V band; R: the B band; I: the I band; O: the optical band; X: the X-ray band; Column 
6 : the references for columns 4 and 5; Column 7: the Kerr black hole masses estimated from 
equation ( 6 a); Column 8 : the Schwarzchild black hole masses estimated from equation ( 6 b); 
Column 9: the black hole masses estimated with other methods in the literatures; Column 10: 
the references for column 9. 
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